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ABSTRACT 

Estimates  are  made  of  the  energy  required  to  expand  an  exploding  shell 
casing  relative  to  the  kinetic  energies  of  the  fragments.  It  is  shown  that 
the  energy  of  deformation  detracts  very  little  from  the  ultimate  fragment 
velocities. 
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ATTENUATION  OF  FRAGMENT  VELOCITIES  DUE  TO  THE  DEFORMATION  RESISTANCE 

OF  THE  CASING 

In  predicting  the  velocity  of  fragments  produced  by  an  exploding  shell, 
neither  Gurney1  nor  Thomas2  take  into  consideration  the  energy  required  to 
expand  the  metal  casing  against  the  resulting  tensile  stresses  which  produce 
the  plastic  flow.  Observations  that  the  fragment  velocities  obtained  with 
relatively  brittle  materials,  such  as  cast  iron,  are  essentially  the  same  as 
those  obtained  with  ductile  casing  materials  would  indicate  that  this  energy 
is  small  compared  with  the  kinetic  energies  imparted  to  the  fragments  and  the 
explosion  products. 

The  following  approach  to  an  analytical  evaluation  of  the  restraining 
effect  is  by  no  means  rigorous.  The  model  it  employs  is  that  of  a  gas, 
initially  at  a  very  high  pressure,  expanding  uniformly  within  the  confines  of 
a  thin-walled  cylindrical  casing  of  infinite  length.  It  does  not  consider  any 
of  the  details  of  the  detonation  process.  Although  the  model  is  not  exact,  it 
may  provide  a  basis  for  evaluating  the  orders  of  magnitude  of  the  various 
factors  contributing  to  the  final  fragment  velocities. 

Consider  first  the  conditions  of  velocity  and  pressure  existing  within 
the  explosive  gases.  At  any  time  t  the  velocity  of  the  gases  may  be  assumed  to 
vary  linearly  from  zero  along  the  axis  of  the  casing  Cr=°)  to  the  value  vR 
adjacent  to  the  casing  (r=R) ,  as  given  by 

v  *  £  v  (1) 

Since  both  vR  and  R  depend  upon  t,  the  gas  velocity  will  depend  upon  both 
r  and  t. 

In  order  for  the  gases  to  accelerate  as  they  expand,  a  pressure  gradient 
must  exist  which  results  in  a  net  force  acting  radially  outward  on  each  element 
of  gas  volume,  as  illustrated  in  Figure  1.  Equating  this  net  force  to  the  mass 
per  unit  length  of  the  element  times  its  radial  acceleration  results  in 

/<VRo\  d2r 

[p(r,t)  -  p(r  dr,t)]  rd0  *( - 1  rdrd0  — —  (2) 

\  *R2  /  dt2 


where  p(r,t)  «  pressure  at  radius  r  at  time  t, 

Pc  =  initial  mass  density  of  the  explosive, 
Rq  =  initial  radius  of  the  casing. 


1  GURNlBY,  R.  W.  The  Initial  Velocities  of  Fragments  from  Bombs,  Shell,  and  Grenades.  U.  S.  Army 
Ballistic  Research  Laboratories,  Aberdeen  Proving  Ground,  Report  No.  405,  September  1043. 
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THOMAS,  L.  H.  U.  S.  Army  Ballistic  Research  Laboratories,  Aberdeen  Proving  Ground,  Report  No.  475. 
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Figure  I.  FORCES  ACTING  WITHIN  AN 
EXPLODING  SHELL 


p(r  ♦  dr,t)  -  p(r,t)  ♦  ||  dr. 


Equation  2  becomes 

O  r2  2 

_  ip.  x  ^£_o  d_r  . 
3r  R2  dt2 

This  may  be  rewritten 


3p  /Ro\  (dv  dr  3v  dR\ 
3r  *  pc \jr)  \3r  dt  4  5R  dt)‘ 


(3) 


(4) 


(5) 

(6) 


Prom  Equation  1 


and 


3v  VR 
3r  *  R 

3 v  r  dvR  r 
3R  *  R  dR  ‘  r2  V 


(7a) 

(7b) 


Substituting  in  Equation  6, 


or 


.  iE 

3r 


p  _R_ 
c  o 


R3  dR 


-  -  2  f  (R)r, 


(8a) 

(8b) 


where  f(R)  is  defined  by 


f(R) 


R2  vn  dv 
o  R 

2  R35 


R 


Equation  8b  can  be  integrated  to  give 

p(r,t)  =  p  (o , t)  -  f(R)r2. 


(9) 


(10) 


In  order  to  evaluate  the  term  p(o,t),  the  equation  of  state  for  the  adiabatically 
expanding  gases  is  taken  to  be 


2 


a 


or 


(11a) 

(lib) 


p  VY  =  p  VY 

Kav  *o  o 

p  R2y  =  p  R2^, 

*av  1  o  o 

where  V  *  instantaneous  gas  volume, 

Vq  =  initial  gas  volume, 

pQ  =  pressure  the  instant  after  detonation, 

pav  »  an  average  pressure  within  the  casing,  which  is  taken  to  be  the 
average  of  p(o,t)  and  p(R,t)  or,  from  Equation  10 


pav  ■  p(o.t)  -  1/2  f(R)  R2. 


(12) 


Combining  Equations  11  and  12  and  substituting  for  p(o,t)  in  Equation  10 


Hence,  the  pressure  acting  against  the  casing  is 

p(R.t)  ■  p/-^  -  y  f(R)R2.  (14) 

The  internal  pressure  given  by  Equation  14  produces  a  tensile  stress  in 
the  cylindrical  casing  which  opposes  the  expansive  effect  of  the  pressure. 

The  forces  acting  per  unit  length  on  an  element  of  the  casing  are  shown  in 
Figure  1,  where  o  is  the  tensile  stress  resisting  the  influence  of  the  inter¬ 
nal  pressure.  The  sum  of  the  radial  components  of  these  forces  may  be  set 
equal  to  the  mass  of  the  casing  element  times  its  radial  acceleration. 

Thus , 

h2d 

p(R,t)Rd0  -  ohde  *  pMhRd6  (15) 

where  h  =  instantaneous  thickness  of  the  casing, 

p  =  mass  density  of  the  casing  material. 

M 


Conservation  of  volume  requires  that 

h-(VR)V  (l6) 

where  h0  is  the  initial  thickness  of  the  casing.  Substituting  for  p(R,t) 
from  (14) ,  Equation  15  may  be  written 
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r 

% 


(17) 


R  2V(2Y  -  «)  -  i 


o  o 


,  dvR  Ro  dVR 

4  pcRo  VR  dlT  ‘  oho  IT  “  pMhoRovR  W‘ 


Johnson  et  al3  have  studied  the  stress-strain  relationship  for  several 
metals  during  high  deformation  rates  and  report  basically  two  types  of  behavior: 
(1)  the  stress-strain  curve  is  similar  to  that  found  in  a  conventional  or 
static  tensile  test  but  with  a  rate  of  work  hardening  roughly  twice  the  static 
value;  and  (2)  the  stress  jumps  immediately  to  a  high  value  in  excess  of  the 
static  tensile  strength  and  remains  constant  thereafter,  independent  of  the 
strain.  For  the  first  case  a  relationship  of  the  form 

o  *  °0En  (18) 

may  be  assumed  where 

e  .  In  |-  (19) 

0 

and  o0  and  n  may  be  taken  to  be  twice  their  static  values.  The  second  case 
would  be  described  by 


a  *  U 


(20) 


where  U  may  be  estimated  as  three  times  the  conventional  ultimate  tensile 
strength.  Normally  n  is  expected  to  decrease  with  increasing  strain  rate, 
contrary  to  the  above  report.  However,  the  behavior  described  by  Johnson  et  al 
will  be  considered  since  it  represents  a  more  extreme  situation  and  will  yield 
a  conservatively  high  estimate  of  the  deformation  energy. 


For  the  first  case  (17)  becomes 

dv, 


p0ro21,r  (2Y  " l3  -  t  pcro  vr  ar  •  °o(ln  *0)  ho  r  *  ’AVr  air 


which  may  be  integrated  directly  to  give 

-2e(y  -  1) 


P  R 
*o  o 


.  1  -  e 


2(Y  -  1) 

Similarly  the  second  case  becomes,  after  integration. 


'1  °oho  n  ♦  1  1  /  .  1  _  \  2 

J’lTTT  e  "  T^o  +  T  pcRo)vR* 


2(y 


R 

L_2 _  j 

-  1)  1 


1  .  e"2e(Y  "  15|-  UhQe  -  j 


’]  -  “V  ■  7  *  7  pcRo)vR: 


Multiplying  by  2ttRq  both  Equations  22  and  23  can  be  put  in  the  form 


(21) 


(22) 


(23) 


SJOHNSON,  R.C.,  STEIN,  B.  A.,  and  DAVIES,  R.  S.  Basic  Parameters  of  Metal  Behavior  Under  High  Rate 
Forming.  Arthur  D.  Little,  Inc.,  Contract  DA-19-020-ORD-5239,  Fifth  Interim  Report  to  U.  S.  Army 
Materials  Research  Agency,  WAL  TR  111.2/20-4,  April  1962  -  June  1962. 
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(24) 


*(m  ♦  i  c)v*  -|Pe(T°.  D  [1 

where  M  ■  the  mass  of  the  casing, 

C  *  thenmass  of  the  explosive, 
ooe 

o'  *  — T  for  o  «  a  en 

n  ♦  1  o 

o'  *  U  for  a  ■  U. 

Gurney's  formula  is  identical  in  form: 


where  v  is  the  ultimate  velocity  of  the  casing  and  E  is  the  energy  per  unit 
mass  of  explosive  ultimately  converted  to  the  kinetic  energies  of  the  casing 
and  the  explosive  gases.  Equation  24  may  be  rewritten 


where  E  is  the  energy  per  unit  mass  of  explosive  which  has  actually  been  con* 


verted  to  kinetic  energy  as 

a  function  of  time.  Thus, 

E  ■  _ Eo_ 

PC(Y  - 

(27) 

and  for  long  times 

E  ■  E. 

(28) 

The  first  term  in  Equation  27  is  the  contribution  of  the  internal  gas 
pressure  to  the  kinetic  energies  of  the  casing  and  explosive  gases.  It  increases 
with  strain  (and  hence  time)  and  approaches  a  final  limiting  value  at  large  val¬ 
ues  of  the  strain.  At  a  strain  of  0.7  the  casing  and  gases  should  have  achieved 
over  90  percent  of  their  final  kinetic  energies.  This  is  exactly  the  strain  (or 
time)  dependence  predicted  by  Thomas  (Ref.  2,  Equation  5.7).  In  reality,  however^ 
a  final  limiting  value  is  probably  reached  at  an  earlier  stage  than  anticipated* 
from  (27)  since,  when  the  casing  fractures  (or  shortly  thereafter),  the  effective¬ 
ness  of  the  expanding  gases  will  be  reduced  due  to  the  gases  escaping  around  the 
fragments.  Hence  the  kinetic  energy  actually  derived  from  the  internal  gas  pres¬ 
sure  will  be  less  than  that  computed  from  (27). 

The  second  term  in  (27)  is  a  negative  contribution  to  the  kinetic  energy 
arising  from  the  resistance  of  the  casing  to  plastic  flow.  As  would  be  expected 
this  term  is  larger  for  small-diameter,  thick-walled  casings,  i.e.,  larger  val¬ 
ues  of  h0/R0.  Furthermore,  since 

_£o  .  2  PM/pc,  (29) 

Ro  c/M 
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the  resistance  to  plastic  flow  has  a  greater  effect  when  the  charge: mass  ratio 
is  low.  The  term  increases  continuously  as  the  strain  increases*  implying  that 
if  the  casing  were  capable  of  expanding  indefinitely  without  breaking*  its  kin¬ 
etic  energy  would  eventually  revert  entirely  to  energy  of  deformation  and  the 
expansion  of  the  casing  would  cease.  In  reality*  however*  even  ductile  casing 
materials  will  fracture  before  the  radius  has  expanded  to  1-1/2  times  its 
initial  size  (corresponding  to  a  strain  of  about  0.4). 

The  extent  to  which  the  resistance  of  the  casing  to  plastic  flow  detracts 
from  the  energy  available  from  the  expanding  gases  may  be  seen  in  Table  I  for 
both  types  of  stress-strain  behavior.  Here  are  tabulated,  for  various  values 
of  the  strain  e,  the  two  terms  appearing  in  Equation  27  computed  for  three  dif¬ 
ferent  size  casings  and  two  values  of  the  initial  gas  pressure.  The  casing 
material  is  taken  to  be  AISI  1060  steel*  quenched  and  tempered  at  1200  F*  for 
which  the  following  constants  apply: 

oQ  *  3.7  x  10s  psi* 

n  «  0.3* 

U  *  400,000  psit 
pM  ■  0.28  lb/cu  in. 

Table  I.  CONTRIBUTIONS  TO  THE  KINETIC  ENERGY  OF  AN  EXPLODING  CASING 


€ 

Energy  Contribution 
from  Internal 

Pressure  (in-lb) 

Pq  <Psi> 

Negative  Energy  Contribution  fron  Resistance  to  Plastic  Flow 

(in-lb) 

cr  ■  CToC11 

Ro/ho 

a  ■  U 

RQ/ho 

3  x  10® 

2  x  106 

10 

5 

3 

10 

5 

3 

0.01 

3.7  x  108 

2.5  x  108 

0.01  X  108 

0.02  x  108 

0.03  x  108 

0.01  x  108 

0.03  x  108 

0.04  x  108 

o.os 

17.6 

10.7 

0.08 

0.15 

0.25 

0.07 

0.13 

0.20 

0.1 

32.3 

21.6 

0.18 

0.36 

0.60 

0.13 

0.26 

0.39 

0.2 

55.0 

36.6 

0.45 

0.89 

1.49 

0.26 

0.53 

0.79 

0.3 

71.0 

47.4 

0.76 

1.51 

2.52 

0.39 

0.79 

1.18 

0.04 

82.3 

54.8 

1.10 

2.20 

3.67 

0.53 

1.05 

1.58 

Fracture  —  —  . 

-  _  _  , 

.  —  -  — 

.  -  — _ 

0.5 

(90.2) 

(60.2) 

1.10 

2.20 

3.67 

0.53 

1.05 

1.58 

0.6 

(95.8) 

(63.8) 

1.10 

2.20 

3.67 

0.53 

1.05 

1.58 

0.7 

(99.7) 

(66.4) 

1.10 

2.20 

3.67 

0.53 

1.05 

1.58 

00 

(109.) 

(72.6) 

1.10 

2.20 

3.G7 

0.53 

1.05 

1.58 

*Two  times  static  test  values. 
tThree  times  static  test  value. 
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The  explosive  was  considered  to  be  Composition  B,  for  which  the  following  con¬ 
stants  were  assumed: 

pc  =  0.061  lb/cu  in 
y  =  2.75 

p0  =  3  x  106  psi  (Case  1) 

Pq  =  2  x  106  psi  (Case  2). 

It  must  be  realized  that  the  negative  contribution  to  the  kinetic  energy 
due  to  the  resistance  to  plastic  flow  can  get  no  larger  than  those  values  which 
correspond  to  the  strain  at  the  time  of  fracture.  Hence,  the  maximum  energy 
absorbed  would  be  that  associated  with  the  maximum  attainable  strain,  i.e., 
e  *  0.4.  Table  I  is  representative  of  casings  exhibiting  this  ductile-type 
behavior.  Brittle  casings  would  fracture  at  a  much  lower  value  of  strain  and 
would  therefore  absorb  considerably  less  energy  than  ductile  casings.  Thus, 
from  Table  I,  a  brittle  casing  which  might  fracture  at  a  strain  of  0.1  would 
absorb  less  than  25  percent  of  the  energy  a  ductile  casing  would  absorb. 

Even  for  the  most  unfavorable  conditions,  namely  a  thick-walled  ductile 
casing  (R0/hQ  =  3);  the  lowest  initial  pressure  considered  to  be  reasonable 
(p0  =  2  x  10*>  psi);  and  the  stress-strain  behavior  described  by  Equation  18, 
the  detrimental  influence  of  the  energy  of  deformation  amounts  to  less  than 
7  percent  of  the  energy  contributed  by  the  internal  gas  pressure  up  to  the 
instant  of  fracture.  In  reality,  the  gases  may  continue  to  effectively  accel¬ 
erate  the  fragments  even  after  fracture  has  occurred,  making  the  values  given 
in  Table  I  for  the  energy  contributed  by  the  internal  pressure  valid  beyond  the 
point  of  fracture.  Thus,  the  energy  of  deformation  would  negate  an  even  smaller 
percentage  than  the  previously  cited  energy  derived  from  the  gas  pressure. 
Furthermore,  a  7  percent  change  in  the  energy  available  produces  just  slightly 
more  than  2.5  percent  change  in  the  fragment  velocities.  Such  a  small  difference 
could  easily  go  undetected  or  could  be  attributed  to  other  variables  in  an  actual 
fragmentation  test. 

The  fragment  velocity  is  plotted  in  Figure  2  as  a  function  of  radius,  strain, 
or  time  (all  are  related)  for  a  charge :mass  ratio  equal  to  one  (corresponding  to 
R0/h0  equal  to  about  10)  and  an  initial  pressure  of  2  x  106  psi.  The  ultimate 
velocity  agrees  favorably  with  that  calculated  from  Gurney's  formula,  where 
E  *  55.7  x  108  inch-pounds  for  Composition  B.  Better  agreement  could  be  obtained 
by  suitable  adjustments  in  p0  and/or  y,  but  such  manipulations  are  probably  not 
meaningful  in  view  of  the  approximations  and  assumptions  already  employed. 

The  conclusion  to  be  drawn  is  that  the  resistance  of  the  casing  to  plastic 
flow  causes  only  minor  attenuations  in  the  fragment  velocities.  The  energy 
attenuation  is  inversely  proportional  to  the  ratio  of  the  shell  radius  to  wall 
thickness,  and  amounts  to  a  reduction  of  about  2.5  percent  in  the  fragment 
velocities  when  the  value  of  this  ratio  is  as  low  as  three. 
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It  is  interesting  to  note  the  difference  in  fragment  velocities  at  the 
time  of  fracture  for  ductile  and  brittle  casings.  The  ductile  casing  is  already 
quite  close  to  its  ultimate  velocity  when  fracture  occurs,  whereas  the  fragments 
from  the  brittle  casing  must  increase  their  velocities  by  nearly  75  percent  in 
order  to  achieve  the  same  velocity.  This  disparity  could  be  overcome  if  the  ex¬ 
plosive  gases  were  to  continue  to  accelerate  the  fragments  after  fracture  in 
spite  of  a  tendency  for  the  gases  to  escape  through  the  spaces  between  the 
fragments.  The  requisite  forces  could  be  derived  both  from  a  normal  component 
due  to  the  retained  internal  pressure  and  from  a  shear  component  resulting  from 
the  drag  forces  of  the  escaping  gas,  Gurney4  has  made  some  quantitative  calcu¬ 
lations  which  support  this  point  of  view. 


*GURNEY,  R.  W.  U.  S.  Army  Ballistic  Research  Laboratories,  Aberdeen 


Proving  Ground,  Report  No.  635. 
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